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t-H ' Abstract. Aims. Non-local thermodynamical equilibrium (NLTE) line formation for neutral and singly-ionized calcium is 

considered through a range of spectral types when the Ca abundance varies from the solar value down to [Ca/H] = -5. We 
evaluate the influence of departures from LTE on Ca abundance determinations and inspect the possibility of using Ca i / Ca n 
line-strength ratios as indicators of surface gravity for extremely metal-poor stars. 

Methods. A comprehensive model atom for Ca i and Ca n is presented. Accurate radiative and electron collisional atomic data are 
incorporated. The role of inelastic collisions with hydrogen atoms in the statistical equilibrium of Ca i/ii is estimated empirically 
from inspection of their different influence on the Cai and Can lines in selected stars with well determined stellar parameters 
^ ■ and high-quality observed spectra. 

Results. The dependence of NLTE effects on the atmospheric parameters is discussed. Departures from LTE significantly affect 
the profiles of Ca i lines over the whole range of stellar parameters considered. However, at [Ca/H] > -2, NLTE abundance 
correction of individual lines may be small in absolute value due to the different influence of NLTE effects on line wings and the 
line core. At lower Ca abundances, NLTE leads to systematically depleted total absorption in the line and positive abundance 
corrections, exceeding +0.5 dex for Ca i A 4226 at [Ca/H] = -4.9. In contrast, NLTE effects strengthen the Ca n lines and lead to 
negative abundance corrections. NLTE corrections are small, < 0.02 dex, for the Can resonance lines over the whole range of 
stellar parameters considered. For the IR lines of multiplet 3d - 4p, they grow in absolute value with decreasing Ca abundance 
exceeding 0.4 dex in metal-poor stars with [Fe/H] < -3. As a test and first application of the Ca i/ii model atom, Ca abundances 
are determined on the basis of plane-parallel LTE model atmospheres for the Sun, Procyon (F IV-V), and seven metal-poor 
, stars, using high S/N and high-resolution spectra at visual and near-IR wavelengths. Lines of Cai and Can give consistent 

abundances for all objects (except Procyon) when collisions with hydrogen atoms are taken into account. The derived absolute 
K> ' solar Ca abundance (from Ca i and Ca n lines) is log £ Ca , = 6.38 ± 0.06. For Procyon, the mean Ca abundance from Ca i lines is 

markedly subsolar, [Ca/H] = -0. 14+0.03. All metal-poor stars within our sample show an overabundance of calcium relative to 
iron with [Ca/Fe] abundance ratios of 0.26 to 0.46 typical for the halo population. The W(Cai4226) / VK(Ca n8498) equivalent 
width ratio is predicted to be sensitive to surface gravity for extremely metal-poor stars, while this is not the case for the ratio 
involving the Can resonance line(s). 
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1. Introduction 

Calcium is one of the best observable chemical elements in 
late-type stars. The resonance lines of Ca i at A 4226 A and Ca n 
at A 3933 and 3968 A lie in the visual spectral range and can 
be measured even in extremely metal-poor stars with metallic- 
ity [Fe/H] < -5 (Christlieb et al. 2001 Frebel et al. l2005l In 
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such stars, Ca is the only chemical element which is visible in 
two ionization stages, and the Ca i and Ca n lines can be potent 
tools in the derivation of accurate values for fundamental stel- 
lar parameters and for the Ca abundance itself. Calcium is an 
important chemical element to study the history of a-process 
nucleosynthesis in the Galaxy. The subordinate lines of neu- 
tral Ca located in the relatively uncrowded yellow-to-red spec- 
tral regions are suitable for spectroscopic analysis over a wide 
range of Ca abundance from super-solar values down to [Ca/H] 
= -4 (Cayrel et al. I2004t . The subordinate lines of ionized 
Ca at A 8498, 8542, and 8662 A are among the strongest fea- 
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tures in the near-infrared spectra of late-type stars with metal- 
licity down to [Fe/H] = -3 (Mallik \TW]\ . The Can triplet 
lies at the focus of cool star research with the advent of large 
spectroscopic surveys of the Galaxy like the Radial Velocity 
Experiment (RAVE, Steinmetz et al. 2006 1 and the upcoming 
ESA Gaia satellite mission (Perryman et al. 2001). These broad 
lines are also powerful abundance and metallicity indicators out 
to large distances via medium-resolution spectroscopy, favor- 
ably located near the flux maximum of red giants. Applications 
range from quantitative analyses of individual stars in glob- 
ular clusters (e.g. Armandroff & da Costa [T9"9"T1 Rutledge et 
al. [T997 1 to stars and stellar populations in nearby galaxies 
(e.g. Tolstoy et al. 1200 II Ibata et al. 12005b . Moreover, the Can 
triplet is useful for the analysis of unresolved stellar systems 
like early-type galaxies (e.g. Saglia et al. 2002). However, as 
an increasing number of high-quality (echelle) spectra of stars 
and unresolved stellar systems becomes available, an equally 
high level in the theoretical modelling is required. 

Previous investigations of the statistical equilibrium (SE) 
of neutral Ca in the Sun and Procyon (Watanabe & Steenbock 
119851 and in the models with effective temperatures T e s = 
4500 K - 6200 K, surface gravities log g = 4.5 -1.0 and metal 
abundance, [Fe/H] = and -1 (Drake [T991l) found rather small 
departures from LTE. Drake J 1 99 li t concluded that "overion- 
ization effects (for Ca i) are not expected to increase in severity 
in metal-poor stars". NLTE analysis of neutral Ca in the sam- 
ple of 252 dwarf and subgiant stars was performed by Idiart & 
Thevenin (2000 1. Previous NLTE investigations for Can were 
concerned only with the resonance lines in the Sun (we cite 
only the first paper, Shine & Linsky [1974b and the IR lines 
of multiplet 3d - Ap (A 8498, 8542 and 8662 A) in moderately 
metal -poor stars by J0rgensen et al. ( 1992 1 for [Fe/H] > — 1 and 
by Andretta et al. (2005) for [Fe/H] > -2. 

The need for a new analysis is motivated by two points. 
Firstly, a fairly extensive set of accurate atomic data on pho- 
toionization cross-sections and oscillator strengths was re- 
cently calculated in the Opacity Project (OP; see Seaton et al. 
( 1994 1 for a general review). We include in the model atom 66 
terms of Cai, contrary to 16 terms in the works of Watanabe 
& Steenbock ( 11985b and Drake ( 11991b . and 36 terms of Can, 
contrary to 3 terms considered by J0rgensen et al. ( 1992) and 
Andretta et al. ( 120051 ). This allows to investigate the detailed 
line formation of the high excitation lines of Cai and Can. 
All our results are based on line-profile analysis and take ad- 
vantage of the advanced theory of collisional broadening by 
atomic hydrogen treated recently by Anstee & O'Mara ( 1995 ), 
Barklem & O'Mara Jl997l [T99%b and Barklem et al. < fl998l . 
Hereafter, these four important papers are referred to collec- 
tively as A&O'M. The second point is connected with an exten- 
sion of the metallicity range observed in stars down to [Fe/H] 
= -5.45 (Aoki et al. 2006). It is highly desirable to have an un- 
derstanding of the likely influence of any departures from LTE 
over this whole range. In stars differing significantly in their 
metal abundance, different subsets of Ca lines are used to cal- 
culate the Can/Cai ionization equilibrium and Ca abundance. 
We investigate here NLTE line formation of our extended list 
of Ca i and Ca n lines in the metallicity range between [Fe/H] 
= and -4.34. An exception is the resonance lines of Cai and 



Ca ii and the Ca n lines of multiplet 3d - Ap in solar metallicity 
stars: their cores and inner wings are most probably influenced 
by the chromospheric temperature rise and a non-thermal and 
depth-dependent chromospheric velocity field neither of which 
are part of the homogeneous photospheric models used in this 
study. 

For the statistical equilibrium of atoms in cool stars, an 
important issue debated for decades, from Gehren d 1975b to 
Belyaev & Barklem ( 2003 ), is the role of inelastic collisions 
with hydrogen atoms. One of our aims is to empirically con- 
strain the efficiency of this type of collisions in the SE of Ca i/ii 
from inspection of their different influence on the Cai and 
Can lines in the Sun, Procyon and seven metal-poor (-1.35 > 
[Fe/H] > -2.43) stars with well determined stellar parameters 
and high quality observed spectra. 

The paper is organized as follows. The model atom of 
Cai/n and atomic data are described in Sect.|2] Departures 
from LTE for Ca i and Ca n, NLTE trends with metallicity, ef- 
fective temperature and surface gravity, and errors of NLTE 
results caused by the uncertainties of atomic data are dis- 
cussed in Sect.|3 In Sect.@J the solar Ca spectrum is studied 
to provide the basis for further differential analysis of stellar 
spectra. We determine the absolute solar Ca abundance using 
four different sets of oscillator strengths based on laboratory 
measurements, OP calculations, and collected from the NIST 
(http://physics.nist.gov/PhysRefData) and VALD (Kupkaet al. 
1999) databases and, thus, estimate the accuracy of the avail- 
able atomic data. Observations and stellar parameters of our 
sample of stars are described in Sect.[5] For the selected stars, 
Sect.[6]investigates whether or not Ca abundances derived from 
the Ca i and Ca n lines agree. NLTE formation of the Ca i res- 
onance line in the atmospheres of metal-poor stars is tested in 
Sect.0 The predictions are made in Sect.[8] for the Ca i/Ca n 
line-strength ratios as indicators of surface gravity for ex- 
tremely metal-poor stars. Our recommendations and conclu- 
sions are given in Sect.[9] 

2. The method of NLTE calculations for Cai/n 
2.1. The model atom 

Energy levels. Calcium is almost completely ionized through- 
out the atmosphere of stars with effective temperature between 
5000 K and 6500 K, with a small fraction (no more than several 
parts in a thousand) of Cai in line formation layers. Minority 
species are particularly sensitive to NLTE effects because any 
small change of the ionization rates largely changes their pop- 
ulations. In order to provide close collisional coupling of Cai 
to the continuum electron reservoir and consequently establish 
a realistic ionization balance between the atomic and singly- 
ionized species, the atomic model for calcium has to be fairly 
complete. Energy levels up to 0.17/0.67 eV below the ioniza- 
tion threshold are included explicitly in our Ca i/n model atom. 
Only the ground state of Cam is considered. The Cai levels 
belong to singlet and triplet terms of the Ami (n < 9 and 
/ < 4), 3dnl (nl = Ap and 3d), and Ap 2 electronic configura- 
tions. Triplet fine structure is neglected except for the Ap 3 P° 
and 3d 3 D splitting. Singlet and triplet ng, 8/, 9d, and 9f lev- 
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Fig. 1. The Ca i/ii model atom. The Ca i and Ca n spectral lines 
used in Ca abundance analysis arise in the transitions shown as 
continuous lines. 

els are combined into a single level due to their small energy 
differences. The final model atom includes 63 levels of Cai. 
The Ca n levels belong to doublet terms of the nl (n — A - 9 
and I < 5), 10s, and lOp electronic configurations. Fine struc- 
ture splitting sub-levels are included explicitly for the terms 
3d 2 D, 4p 2 P°, and Ad 2 D. ng and nh (n = 7, 8, and 9) levels are 
combined into a single level. The final model atom contains 37 
levels of Ca n. 

The energy levels were taken from the NIST atomic spec- 
tra database (Sugar & Corliss 1985 1 for all the Cai electronic 
configurations with n < 5 and Can electronic configurations 
with n < 6. For the remainder, Opacity-Project data from the 
TOPBASE database (Cunto & Mendoza[T592i are used. The 
corresponding term diagram is shown in Fig.^ 

Radiative transitions. Our NLTE calculations rely heavily 
on the OP radiative data. They cover the whole range of al- 
lowed transitions in our model atom. Currently available OP 
calculations do not include the effects of spin-orbit interaction, 
and only average-multiplet oscillator strengths are reported. 
Where required, we decompose the LS multiplet averages into 
their LSJ fine structure components in correspondence with 
their relative line strengths. Additional radiative coupling be- 
tween the different spin systems of Ca i is provided by the tran- 
sition As 2 - 4p 3 P° (Drozdowski et al. 119971 and by all inter- 
combination transitions listed by Smith & Raggett ( 1981 ). The 
radiative transfer equation is solved, in total, for 421 line tran- 
sitions in Cai and 213 transitions in Can. 

We compared OP fj for Ca i with the experimental values 
from S habano va {T963>, Smith & Gallagher < fl966b , Smith & 
O'Neil ( fT9"73l . Smith & Raggett fiWTl and Smith fT9"S8l . in 
total for 24 multiplets with log gf ranging between 0.73 and 



-2.25. The difference A log gf (laboratory - OP) shows no cor- 
relation with the multiplet strength, and the mean value equals 
-0.07 + 0.28. For Ca n, measured oscillator strengths are avail- 
able only for the multiplets As - Ap and 3d - Ap (Theodosiou 
1989). They agree within 12% with the corresponding OP val- 
ues. OP data for the remaining Ca i and Ca n transitions were 
compared with calculations of Kurucz J1992i . For the 151 
strongest (logg/(OP) > -2.5) Cai multiplets, the mean dif- 
ference A log gf (Kurucz - OP) = -0.03+0.93. Kurucz calcu- 
lations predict, on average, lower transition probabilities com- 
pared to OP data for the 58 strongest (log gf (OP) > -2.5) Can 
multiplets, and the mean difference A log gf (Kurucz - OP) = 
-0.17±0.86. With respect to systematical errors OP data for 
Cai and Can transitions are expected to be accurate within 
15%. Additional arguments for the validity of this statement 
come from the analysis of solar Ca spectrum (see Sect.@}. 

Photoionization from all energy levels is treated utilizing 
OP cross-sections as available through the TOPBASE database. 
Cross-sections for the combined ngh levels in Ca n are assumed 
to be equal to that for the corresponding ng levels. 

Collisional transitions. In our calculations we take into ac- 
count inelastic collisions with electrons and hydrogen atoms 
leading to both excitation and ionization. Drawin's ( 1968) for- 
mula as described by Steenbock & Holweger ( 1984) is widely 
used to calculate hydrogenic collisions, and it suggests that 
their influence is comparable to electron impact. Recently it 
was shown both experimentally (Belyaev et al. 1999 for Na I) 
and theoretically (Belyaev & Barklem 2003 for Li I) that 
Drawin's formula overestimates the collision rate of the res- 
onance transitions by several orders of magnitude. However, 
for SE calculations, estimates are required for transitions be- 
tween all states which might affect the population of the states 
of interest. In this study, we constrain the efficiency of hydro- 
genic collisions empirically. It is represented by a scaling fac- 
tor Sh applied to Steenbock & Holweger's formula. The cross- 
sections calculated using this formula were multiplied by Sh — 
(no hydrogenic collisions), 0.01, 0.1, and 1 in order to agree 
Ca abundances derived from the Ca i and Ca n lines in the Sun 
and selected stars. 

For electron impact excitation, detailed results are avail- 
able from the /^-matrix calculations of Samson & Berrington 
(2001) for the transitions from the ground state to all the ex- 
cited levels up to 3d4p P° in Cai and from calculations of 
Burgess et al. IT995) in a non-exchange distorted wave ap- 
proximation for the transitions between the lowest seven terms 
(up to 4/ 2 F°) in Can. Uncertainties in data obtained by the 
/^-matrix method are typically on the order of a few 10%, and 
Samson & Berrington note excellent agreement of their results 
with experimentally derived rates for As 2 -AsAp 'P°. For the re- 
maining bulk of the transitions, approximate formulae must be 
used. Basically, we apply the impact parameter method (IPM) 
of Seaton ( 1962a) to the allowed transitions and the Allen's for- 
mula ( 1973 ) with a collision strength of 1 .0 to the optically for- 
bidden transitions. Electronic collision rates based on Samson 
& Berrington and Burgess et al. data where they are available, 
IPM data and the Allen's formula for the remaining transitions 
are referred to below as the standard collisional recipe. For test 
purposes (see Sect. l3.3> . we treat the allowed transitions absent 
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in the Samson & Berrington and Burgess et al. papers using the 
van Regemorter's formula ( 1962 1. 

Electron impact ionization cross-sections are calculated ap- 
plying the formula of Seaton ( 1962b I with threshold photoion- 
ization cross-sections from the OP data. 

2.2. Model atmospheres. Programs and methodical 
notes 

Calcium is assumed to be a trace element because its contri- 
bution to the continuous opacity and the reservoir of free elec- 
trons is by, at least, one order of magnitude smaller than that 
from the more abundant elements Mg, Si, and Fe. Thus, we ob- 
tain statistical equilibrium populations for Ca i/ii while keeping 
the atmospheric structure fixed. All calculations are performed 
with plane-parallel, homogeneous, LTE, and blanketed model 
atmospheres computed with the MAFAGS code (Fuhrmann et 
al. lTW7l . 

We use a revised version of the DETAIL program (Butler & 
Giddings 1985 1 based on the accelerated lambda iteration fol- 
lowing the efficient method described by Rybicki & Hummer 
(1991 1992 1 in order to solve the coupled radiative transfer 
and statistical equilibrium equations. All bound-bound (b - b) 
and bound-free (b-f) transitions of Ca i and Ca n are explicitly 
taken into account in the SE calculations. The 15 strongest b-b 
transitions in Cai and all LS J transitions As - 4p and 3d - 4p 
in Can are treated using Voigt profile. Microturbulence is ac- 
counted for by inclusion of an additional term in the Doppler 
width. The van der Waals damping parameters based on the 
Anstee & O'Mara's d!995> theory are taken from the VALD 
database (Kupka et al. 1999 1. The remaining b-b transitions 
are treated using depth-dependent Doppler profiles. 

In addition to the continuous background opacity, the line 
opacity introduced by both Hi and metal lines is taken into 
account by explicitly including it in solving the radiation trans- 
fer. The metal line list has been extracted from Kurucz' ( 1992 1 
compilation and contains about 650 000 atomic and molecular 
lines between 1300 A and 300 000 A. Ca lines are excluded 
from the background. 

The obtained departure coefficients are then used to 
compute the synthetic line profiles via the SIU pro- 
gram (www.usm.uni-muenchen.de/people/reetz/siu.html). In 
this step of calculations, Voigt profile functions are adopted 
and the same microturbulence value V m j c as in DETAIL is ap- 
plied. Oscillator strengths and van der Waals damping con- 
stants of the Ca i and Ca n lines and their accuracy are discussed 
in Sect.gl 

Ca is represented in the solar system matter by several iso- 
topes with the isotope abundance ratio 40 Ca : 42 Ca : 43 Ca : 
44 Ca : 46 Ca : 48 Ca = 96.9 : 0.6 47 : 0.135 : 2.09 : 0.004 : 
0.187 (Anders & Grevesse 19891. In this study, we account 
for isotope structure of the Can A 8498 line 1 with the isotope 
shifts measured by Nortershauser et al. (1998 1. Wavelengths 
and the product of gf and solar fractional isotope abundance e 
are given in Tabled for each isotope component. The oscilla- 

1 Of the three IR triplet lines we look at Can A 8498 only, as a) it is 
the weakest of the three and b) it is least perturbed by Paschen lines. 



Table 1. Atomic data for the isotopic components of Can 
A 8498. The fractional isotope abundances e correspond to the 
solar system matter 



A, A 


isotope 




8498.023 


40 


-1.43 


8498.079 


42 


-3.60 


8498.106 


43 


-4.29 


8498.131 


44 


-3.10 


8498.179 


46 


-5.81 


8498.223 


48 


-4.14 



tor strength is taken from measurements of Theodosiou ( 1989 1, 
log gf{A 8498) = -1.416. Isotope shifts are much smaller 
(<10mA) for the resonance lines in Cai and Can (Lucas et 
al. 2004 1 and we treat them as single lines. No data is available 
for the remaining Ca lines. We neglect the hyperfine structure 
of Ca lines due to very low fractional abundance of the only 
odd isotope 43 Ca (0.135%). 

3. NLTE effects as a function of stellar parameters 

In this section, we investigate NLTE effects for Cai/n in the 
following range of stellar parameters: T e g between 5000K and 
6000K, logg = 3.0 and 4.0, and [Fe/H] between and -3. For 
the models with [Fe/H] = 0, the Ca abundance was assumed 
to follow the global metallicity. Ca enhancement with [Ca/Fe] 
= 0.4 was assumed for the metal-poor models. Statistical equi- 
librium was computed using the standard recipe for electronic 
collisions and no hydrogenic collisions were taken into account 
(Sh = 0). 

3.1. Statistical equilibrium of Cai and NLTE effects for 
spectral lines 

Our calculations indicate a general behaviour of departure co- 
efficients, bj = ri; /m , independent of effective tempera- 
ture, surface gravity and metallicity. Here, nf hTE and n[ ,TE are 
the statistical equilibrium and thermal (Saha-Boltzmann) num- 
ber densities, respectively. Departure coefficients of the impor- 
tant levels of Cai are plotted in Fig. |2] for selected models of 
our grid. It can be seen, all levels of Ca i are underpopulated in 
the atmospheric layers above log rjooo = 0. Overionization is 
caused by superthermal radiation of non-local origin below the 
thresholds of low excitation levels of Ca i. In atmospheres with 
solar or mildly deficient Ca abundance, [Ca/H] > — 1, the most 
important levels are 4p 1 P°, 3d 3 D, and 4p 3 P° with the thresh- 
olds at 3898 A, 3450 A, and 2930 A, correspondingly. At the 
lower metallicity and Ca abundance, depopulation processes in 
Ca i are dominated by enhanced ionization of the ground state 
due to a reduction of the continuous absorption coefficient be- 
low the threshold of this level at 2028A. 

Statistical equilibrium of Ca i in the metallicity range be- 
tween and -1 was discussed in detail by Drake (1991 1. 
Referring to the main continuous opacity sources in the near 
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Fig. 2. Departure coefficients log b as a function of log T5000 for 
selected levels of Ca 1, plotted for selected models of the grid. 
For each term of 3af 3 D and Ap 3 P°, departure coefficients of 
the fine structure splitting sublevels are obtained to be close 
together and only the lowest sublevels, 3d i D l and 4p 3 P^, are 
shown. Stellar parameters are quoted in each panel 



ultraviolet wavelength regime, due to the FT ions and excited 
neutral H atoms, Drake has shown that overionization effects 
decrease with increasing model temperature, decreasing sur- 
face gravity, and decreasing Ca abundance. Our calculations 
for the stellar parameters overlapping with that of Drake sup- 
port qualitatively his conclusions. However, quantitatively, in 
our work, overionization effects are found to be smaller. The 
explanation lies with our extended model atom that contains 
the energy levels up to 0.17 eV below the ionization thresh- 
old, contrary to 1.06 eV in Drake's paper. Overionization ef- 
fects on Cai are partially cancelled due to collisional cou- 
pling of the Cai high excitation levels to the ground state of 
the majority species, Can, that keeps thermodynamic equilib- 
rium population. Our model atom provides this closer coupling 
and in this respect is more realistic than the one by Drake. 
The smaller departures from LTE for level populations result 
in weaker NLTE effects for spectral lines and smaller differ- 
ences between derived NLTE and LTE Ca abundance, Anlte = 
log £nlte - log slte- We refer to Anlte as the NLTE abundance 
correction. For example, for atmospheric parameters T e g = 
5800K, log g = 4.5, [Fe/H] = 0, Drake gives nearly equal values 
Anlte = 0.11 dex for ,15261 (multiplet 22), ,16455 (multiplet 
19), and A 6166 (multiplet 20). The corresponding values from 
our calculations are 0.04 dex, 0.05 dex, and 0.07 dex. 

Contrary to the models with [Fe/H] > -1, overionization 
effects increase with increasing model temperature and de- 
creasing global metallicity, when ones goes to the lower metal- 
licity and Ca abundance, [Fe/H] and [Ca/H] < - 1 . The expla- 
nation lies with a behavior of the continuous absorption coeffi- 




6438.6 



6438.8 



6439.0 



6439.2 



6439.4 




6169.40 6169.50 6169.60 6169.70 6169.80 

Fig. 3. NLTE (continuous line) and LTE (dotted line) theoret- 
ical profiles of the selected Cai lines for the models with the 
same T e g = 5500K and log g = 4.0 but with different metallic- 
ities. For the A 6 162 and ,16439 lines, [Fe/H] = 0, -2, and -3 
from bottom to top; for ,16169, [Fe/H] = 0, -1, and -2 

cient below 2028A that is, in metal-poor atmospheres, mainly 
due to quasi-molecular hydrogen absorption (Doyle 1968), the 
H~ ion, and the HJ ion. The contribution of the H ion in- 
creases with increasing effective temperature and decreases 
with decreasing metallicity, while the contribution of quasi- 
molecular hydrogen absorption changes in the opposite direc- 
tion. Their competition results in effects seen in Fig.|2] 

NLTE effects on spectral lines are illustrated in Fig. [3] and 
Table|2] Fig.[5]shows NLTE and LTE profiles of three represen- 
tative lines of neutral calcium computed for the models with 
different global metal abundance. These are the low excitation 
line ,16162 which is the strongest among all subordinate lines, 
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A 6439 with the largest oscillator strength, fy = 0.390, but rela- 
tively small van der Waals damping constant, log C(, = -31.58, 
which is lower by 1.28 dex than that of ,16162, and the inter- 
mediate strength line A 6169.5. Tabled presents NLTE abun- 
dance corrections, A NLTE , for 16 Cai lines. Three other Cai 
lines, ,16161, A 6169.0, and A 6122, have Anlte close to those 
of A 6166, A 6169.5, and ,16162, respectively. NLTE and LTE 
line profiles and equivalent widths were calculated with the 
laboratory oscillator strengths and Ce values given in Table|4] 
Micro turbulence value, V m i c = 1 km s~\ was adopted for all the 
models. Departures from LTE for spectral lines can be under- 
stood by considering the depths of formation of various parts 
of the lines and by inspecting the level departure coefficients 
and line source functions. In the metallicity regime, [Fe/H] 
> -1, our results support the conclusions of Drake dl991> : 
for each investigated line, NLTE effects lead to enhanced ab- 
sorption in the line core and depleted absorption in the line 
wings (Fig.©. This can be understood because the line wings 
are formed in deep layers where overionization depopulates all 
Ca i levels, but the line cores are formed at small depths, above 
log T5000 = -2, where the upper levels of the transitions are 
underpopulated to a more extent than are the lower levels due 
to photon losses in the line wings. The most prominent exam- 
ple is a steep decrease of b(Ap x ~P°) above logrsooo = -2 in 
all metal-poor models (Fig.|2}. The line source function drops 
below the Planck function at these depths resulting in the en- 
hanced absorption in the line cores. Combined effect on the 
line strength is that the NLTE abundance correction is small in 
most cases. Its sign and value depends on the contributions of 
the line core and wings to the overall line strength. NLTE cor- 
rections tend toward negative values with increasing T e ff and 
decreasing log g. They are more negative for the lines of multi- 
plet 18, /1 6439, A 6471, /1 6493, and A 6499 (Table©. All these 
trends reflect the behaviour of the van der Waals broadened line 
wings. 



Contrary to the solar metallicity models, the energy levels 
become weakly coupled far inside the metal-poor atmospheres 
with [Fe/H] < -1 due to deficient collisions (Fig.[2J. For each 
model, at the depths where the weak lines are formed, the up- 
per levels are all depleted to a lesser extent relative to their LTE 
populations than are the lower levels. The lines are weaker rel- 
ative to their LTE strengths not only due to the general overi- 
onization but also due to bjbi > 1 resulting in the line source 
function S /„ b u jb\B v > B v and the depleted line absorp- 
tion. Here, b u and b/ are the departure coefficients for the upper 
and lower levels of the transition. For example, this is valid for 
,16169.5 in the model with T eB = 5500 K, logg = 4.0, [Fe/H] 
= -2 and ,16162 and ,16439 in the model 5500 K / 4.0 / -3. 
The corresponding NLTE and LTE profiles are shown in Fig.0 
If the line is strong for a given set of stellar parameters (e.g., 
A 6162 and ,16439 for the models with [Fe/H] = -2), NLTE 
effects on line profiles are similar to that for the models with 
[Fe/H] > -1. 
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Fig. 4. The same as in Fig.[2]for Can 



3.2. Statistical equilibrium of Ca n and NLTE effects on 
spectral lines 

Figure 0] shows departure coefficients for some important lev- 
els of Can plotted for selected models of our grid. In the tem- 
perature regime we are concerned with here (T e g = 5000 K - 
6000 K), Ca ii dominates the element number density over at- 
mospheric depths. Thus, no process seems to affect the Can 
ground state population, and As keeps its thermodynamic equi- 
librium value. An exception is the uppermost layers above 
logT 5 ooo = -4 in the very metal -poor ([Fe/H] < -2) models 
with T e ff = 6000 K. Can competes there with Cam in contribu- 
tion to element population, and enhanced ionization of the low 
excitation levels of Ca n leads to underpopulation of total Ca n. 
This effect is amplified with decreasing metallicity. The levels 
3d and Ap follow the ground state in deep layers, and their cou- 
pling is lost at the depths where, for each model, photon losses 
in the weakest line A 8498 of the multiplet 3d - Ap start to be- 
come important. These are the uppermost layers above log rsooo 
= -5 in the solar metallicity models and above log rsooo = -3 
in the models with [Fe/H] = -2. At [Fe/H] = -3, detailed bal- 
ance in the transition 3d - Ap is destroyed in the deeper lay- 
ers around log rsooo = -2. The departure coefficients of Ad 2 T) 
and the higher excitation levels begin to deviate from 1 far in- 
side the atmosphere even in the solar metallicity models due to 
photon losses in the transitions to low excitation levels. 

NLTE leads to strengthened Ca n lines and negative NLTE 
abundance corrections. In the following, we consider a be- 
haviour of individual lines. 

In the stellar parameter range we are concerned with in 
this paper, departures from LTE occur only in the very core 
of Can K A 3933. The NLTE abundance correction does not 
exceed 0.02 dex in absolute value. This can be understood be- 
cause the line remains to be strong even at [Ca/H] = -3, and its 
line strength is dominated by the line wings. 
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Table 2. NLTE abundance corrections (dex) for the Ca i and Ca n lines depending on effective temperature, surface gravity, and 
metallicity. Everywhere V m i c = 1 kms -1 . No hydrogenic collisions were taken into account (Su = 0) in SE calculations. Blanks 
denote theoretical NLTE equivalent widths below 5 mA 
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0.12 


0.14 


0.16 


-0.06 




-0.08 


5000/4/-3 


0.23 


0.25 








0.27 




0.12 




0.16 






0.22 




0.24 




-0.16 






5500/3/0 


0.10 


0.05 


0.07 


0.00 


-0.03 


-0.02 


0.06 


-0.07 


0.02 


-0.02 


0.10 


-0.02 


-0.11 


-0.05 


-0.06 


0.00 


-0.02 


-0.20 


-0.22 


5500/3/-2 


0.20 


0.19 


0.32 


11.27 


0.12 


0.10 


0.25 


0.13 


0.1 1 


-0.04 


0.28 


0.23 


-0.02 


0.14 


0.19 


0.18 


-0.11 


-0.02 


-0.16 


5500/3/-3 


0.12 


0.29 








0.29 




0.24 




0.18 






0.24 




0.26 




-0.29 






5500/4/0 


0.08 


0.05 


0.05 


0.01 


0.01 


0.03 


0.05 


-0.01 


0.05 


0.03 


0.08 


0.02 


-0.02 


-0.04 


-0.02 


0.00 


-0.02 


-0.10 


-0.13 


5500/4/- 1 


0.13 


0.08 


0.14 


0.08 


0.06 


0.00 


0.09 


0.00 


0.1(1 


0.00 


0.15 


0.06 


-0.10 


-0.01 


-0.04 


0.04 


-0.03 


-0.09 


-0.21 


5500/4/-2 


0.19 


0.18 


0.28 


0.25 


0.11 


0.10 


0.22 


0.13 


0.1 1 


0.03 


0.25 


0.21 


0.00 


0.12 


0.16 


0.15 


-0.08 




-0.10 


5500/4/-3 


0.20 


(1.27 








0.29 




0.19 




0.18 






0.23 




0.24 




-0.21 






6000/3/0 


0.09 


-0.01 


0.05 


-0.01 


-0.04 


-0.14 


0.02 


-0.13 


0.01 


-0.17 


0.08 


-0.05 


-0.27 


-0.07 


-0.14 


-0.02 


-0.03 


-0.28 


-0.29 


6000/3/-2 


0.10 


0.20 


0.26 


0.21 


0.13 


0.13 


(1.22 


0.16 




-0.01 




0.21 


0.01 


0.1 1 


0.17 


0.14 


-0.16 


-0.04 


-0.20 


6000/3/-3 


0.02 


0.26 








0.27 








(1.19 






0.22 








-0.41 




0.00 


6000/4/0 


0.07 


0.02 


0.04 


0.00 


0.01 


-0.05 


0.01 


-0.05 


0.115 


-0.05 


0.06 


-0.01 


-0.13 


-0.06 


-0.10 


-0.02 


-0.02 


-0.17 


-0.27 


6000/4/- 1 


0.11 


0.03 


0.13 


0.08 


0.05 


-0.10 


0.07 


-0.02 


0.07 


-0.10 


0. 1 3 


0.05 


-0.22 


-0.01 


-0.08 


0.03 


-0.05 


-0.13 


-0.29 


6000/4/-2 


0.13 


0.18 


0.24 


0.20 


0.12 


0.12 


0.19 


0.15 




0.03 




0.19 


0.01 


0.10 


0.14 


0.13 


-0.11 


-0.01 


-0.12 


6000/4/-3 


0.10 


0.25 








0.29 








0.19 






0.21 








-0.23 







Similar to the resonance line, the IR lines of multiplet 
3d - 4p, A 8498, A 8542, and A 8662, reveal NLTE effects only 
in the Doppler core. The line core is strengthened because 
the line source function drops below the Planck function due 
to b u /bi < 1 and the lower level of the transition is over- 
populated (bi > 1) at the line core formation depths. The 
NLTE correction is larger in absolute value for the weakest 
line, A 8498, compared to the other two due to smaller contri- 
bution of the line wings to the overall line strength. However, 
even for A 8498, Anlte is small for the models with [Fe/H] = 
and -1, |A NLTE | < 0.05 (TableEJi. The van der Waals broad- 
ened wings are weakened with decreasing Ca abundance and, at 
fixed metallicity, with decreasing surface gravity and increas- 
ing temperature. NLTE effects are amplified in the same direc- 
tions such that Anlte 8498) reaches -0.41 dex for the model 
6000 K / 3.0 / -3. Significant NLTE effects can be expected to 
occur in the low-density atmospheres of metal-poor red giants. 
Further investigations will be conducted when we have proper 
observations at our disposal. 

The IR high excitation lines of multiplets 5p - 5d, A 8248 
and A 8254, and 4d - 4/, A 8912 and A 8927, can be good can- 
didates for determination of Ca n abundance in close to solar 
metallicity cool stars. They are of intermediate strength and 
nearly free of blends. It can be seen from Fig.0]that the high 
excitation levels are all depleted relative to their LTE popula- 
tions at line formation depths in the models with [Fe/H] > -2. 
However, the lower levels of the considered transitions, 5p and 
Ad, are underpopulated to a lesser extent than are the upper 
levels, 5d and 4/, resulting in smaller line source functions 
compared to the Planck function and enhanced line absorption. 
NLTE effects are significant even for the solar metallicity mod- 
els and become stronger with increasing T e g and decreasing 
logg. The lines of multiplet 5p - 5d are weaker compared to 
A 8912 and A 8927 and, in general, reveal smaller NLTE effects. 
In Table|2 we present Anlte for A 8248 and A 8927. NLTE cor- 
rections for the second line of multiplet 4-d - 4/ are close to 



the corresponding values for A 8927 and, in general, smaller by 
up to 0.03 dex in absolute value. For A 8254, NLTE effects are 
very small with Anlte at the level of a few parts in a hundred. 

3.3. Error estimates for the Ca i/ii SE calculations 

To assess the effects of crucial atomic data on the accuracy 
of NLTE Ca abundances derived from the Ca i and Ca n lines, 
test calculations were performed for three models with T e g = 
5500 K, logg = 4.0 and [Fe/H] = -1, -2, and -3. For each 
parameter or set of cross-sections varied, we computed for a 
given line a small grid at different abundances to determine the 
systematic shift in Ca abundance needed to fit the NLTE equiv- 
alent width evaluated using our standard set of atomic data. The 
results of the tests are summarized for selected lines in Table[3] 

OP photoionization data based on the R-matrix calculations 
typically are accurate to 10%. In test computations, we assume 
a factor of two uncertainty of photoionization cross-sections as 
a worst case. As expected, a variation in the photoionization 
rates affects the lines of the minority species (Ca i) by way of a 
displaced ionization balance and does not affect the lines of the 
majority species (Can). Corrections for the Cai lines decrease 
with decreasing metallicity. The explanation lies with an in- 
creasing fraction of Ca m at the same optical depth in the mod- 
els with decreasing metallicity. When cross-sections of both the 
Ca i and Ca n levels change, the Ca i/Ca n/Ca in ionization bal- 
ance is established in such way that the Ca i/Ca n ratio is influ- 
enced to a lesser extent in the more metal-poor models. 

In the present work, detailed electronic collision excitation 
data are used for a considerable number of transitions in Ca n 
and for all important transitions from the ground state in Ca i. 
In test calculations, we vary collisional rates for the remaining 
transitions applying the van Regemorter's formula (1962) in- 
stead of IPM data. Van Regemorter's collisional rates C,j(Reg) 
are, in general, larger than the corresponding IPM-based val- 
ues, by up to 2 orders of magnitude for the transitions with 
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Table 3. Uncertainties in the NLTE analysis of Ca i/ii. Blanks denote NLTE theoretical equivalent widths below 5 mA 
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energy separation AE < 2 eV (see Fig. 2 in Mashonkina ( 1996) 
for the transitions in Mg i). As a result, NLTE effects are weak- 
ened compared to the standard collisional recipe and a NLTE 
Ca abundance is obtained to be closer to the LTE one (lower 
from the Cai lines and higher from the Can lines, Tabled}- 
The only exception is Cai 5857 for which NLTE effects are 
strengthened. In any case, a variation of electronic collision 
rates only weakly affects the Ca i lines, such that the derived Ca 
abundance changes by 0.04 dex, at maximum. In contrast, the 
effect is significant for the Ca n lines, and a change in Ca abun- 
dance may consist of 0.07 dex to 0.14 dex for different lines. A 
different reaction of the Ca i and Ca n lines can be understood 
because the main mechanism of departures from LTE is con- 
nected with b—f transitions for Ca i, while with b—b transitions 
for Ca n. 

We inspect also the effect of including hydrogenic colli- 
sions in our SE calculations. Table[3] shows the difference in 
Ca abundance derived assuming Su = 1 and ignoring H colli- 
sions (Sh = 0). As expected, departures from LTE are weak- 
ened for the Ca n lines. Changes in Ca n abundance are compa- 
rable to those obtained when varying electronic collision rates. 
A somewhat unexpected behaviour is seen for most Ca i lines 
in the model with [Fe/H] = - 1 and for A 6439 in the model with 
[Fe/H] = -2: NLTE effects seem to be amplified when the to- 
tal collisional rates increase. This can be understood because 
departures from LTE go down in the line wings, as expected, 
but are hardly changed in the cores of strong lines due to the 
inefficiency of collisions in the uppermost atmospheric layers. 
Thus, the line wings do not act anymore to reduce the com- 
bined NLTE effect on the overall line strength, and it becomes 
larger compared to the case of electronic collisions only. 

In addition, we examine the effect of a shift of the line 
formation depths due to increasing the microturbulence value 
in the model. An increase of V m j c by 0.5kms _1 has a neg- 
ligible effect on Anlte for every Ca line in the models with 
[Fe/H] = -2 and -3 and leads to strengthening NLTE effects for 



the model with [Fe/H] = — 1 . The maximum correction amounts 
to -0.03 dex. 

Thus, the largest uncertainty of NLTE results for Ca i/ii is 
caused by poor knowledge of collision processes. Below we 
empirically constrain collisional data by way of analysis of Ca 
lines in solar (SectS and stellar (Sect.|6j spectra. 

4. Solar Ca abundance and Cai /Can ionization 
equilibrium 

In this section, we derive solar Ca abundance from the Cai 
subordinate lines and the Can high excitation lines and ex- 
amine atomic data used in SE calculations and element abun- 
dance determinations. Taking advantage of the applied NLTE 
approach, both weak and strong Ca lines are included in anal- 
ysis. An exception is the resonance lines in Cai and Can and 
the Can lines of multiplet 3d - 4p, because their cores and in- 
ner wings are influenced by the chromospheric temperature rise 
and non-thermal and depth-dependent chromospheric velocity 
field that is not part of the MAFAGS model of the solar atmo- 
sphere. However, we check here the wings of Cai /14226 and 
Can ,18498. 

We use solar flux observations taken from the Kitt Peak 
Solar Atlas (Kurucz et al.[T984i and select the Ca lines free of 
blends, or with only one line wing distorted (e.g., Cai A 4425), 
or with the blending lines which can be correctly taken into 
account (e.g., Cai A6512 in the wing of Balmer line Ha). The 
investigated lines are listed in Table|4] 

Four different sets of oscillator strengths are applied and 
compared in this study, (i) fy obtained from laboratory mea- 
surements. They are available for all selected Cai lines and 
shown in Tabled (column LAB) together with their sources, 
(ii) fij based on OP calculations. They are available for all op- 
tically permitted transitions in Ca i and Ca n, (iii) and (iv) the 
data from NIST and VALD databases, respectively. 

We use the radiative widths obtained by Kurucz J1992I 
from radiative lifetimes and accessible via the VALD database. 
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For 17 Cai lines, the C(, values are computed from damp- 
ing parameters given by Smith ( 1981 1 and based on measured 
parameters for broadening by helium. For the remaining lines, 
we adopt values based on either A&O'M or Kurucz cal- 
culations, giving a preference to the first source. The neces- 
sary data are accessible via the VALD database. The van der 
Waals damping parameters based on the perturbation theory 
of A&O'M agree within 0. 1 dex of log Cc with quantum me- 
chanic computations by Spielfiedel et al. (1991 1 for the Cai 
multiplet 3 (4p 3 P° - 5s 3 S) and by Kerkeni et al. d2004l > for the 
Cai resonance line A 4226. For four among five common Cai 
multiplets, the predicted A&O'M parameters lead to stronger 
collisional broadening compared to that from the experimen- 
tal data of Smith ( 1981 ). The difference in log ranges from 
0.19 dex to 0.32 dex. The opposite is the case for Cai A 5261 
line with the experimental value of log larger by 0.25 dex 
compared to the predicted one. For Cai A 4425 and the lines of 
Ca ii multiplets 5p-5d and 4d-4f, the Ce values were obtained 
empirically from the fitting of solar line profiles. We note that 
not only A 4425 but also two other lines of the Ca i multiplet 
4p 3 P° - 4<i 3 D, A 4454 and A 4455, though being blended, cer- 
tainly cannot be fitted with log C(, = -30.23 computed from 
A&O'M data. We find log C 6 = -30.9 for these lines. The best 
fit of /14425 is shown in Fig|5] For the Can multiplet 5p — 5d, 
log C(, = -30.85 was found from the requirement that element 
abundances derived from the weaker line, A 8254 (W^ =18 
mA), and the stronger line, A 8248 (W^ = 65 mA), must be 
equal. Contrary to multiplet 5p - 5d, the lines of the Can mul- 
tiplet 4d - 4f, A 8912 and A 8927, are both sensitive to a vari- 
ation of the van der Waals damping constant. We assume that 
the broadening parameter calculated by Kurucz is underesti- 
mated to the same extent as the corresponding value for mul- 
tiplet 5p - 5d and, thus, obtain logC6 = —31.1 for multiplet 
4d - 4f. The best fit of A 8927 is shown in Fig|5] The atomic 
data are presented in Tabled 

A depth-independent microturbulence of 0.9 km s 1 is 
adopted. Our synthetic flux profiles are convolved with a pro- 
file that combines a rotational broadening of 1.8kms~' and 
broadening by macroturbulence with a radial-tangential profile 
of Vmac = 3 km s _1 to V mac = 4 km s _1 for different lines. 

For each investigated line, the product logg/e^ a was 
obtained from solar line-profile fitting under various line- 
formation assumptions: NLTE Sh = 0, Sh = 0.1, and Sh — 1 
and LTE. Table@] presents the NLTE values logg/e^ derived 
assuming Sh = and abundance corrections Ax = log ex - 
log £s H =o where X takes the meaning of 0. 1, 1, and LTE for Sh 
= 0.1, Sh = 1, and the LTE assumption, respectively. By def- 
inition, Alte (Table 4) = -Anlte (Table 2). It is interesting to 
note that, for every investigated line, Anlte given in the first 
string of Table[2]does not coincide in absolute value with Alte 
in Table|4] This can be understood because NLTE corrections 
in Table[2]were calculated from comparison of NLTE and LTE 
equivalent widths, while abundance corrections in Tabled are 
based on the analysis of line profiles. For each spectral line, 
the best NLTE fit, independent of Sh value, reproduces the line 
core better than the LTE one. However, for many lines, even 
the NLTE fit is not perfect. As an example, we show in Fig.|5| 
the best NLTE (Sh - 0) and LTE fits for two lines of Cai and 
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Fig. 5. Synthetic NLTE (Sh = 0, continuous line) and LTE (dot- 
ted line) flux profiles of the Cai ,14425, ,16122, and ,16439 
and Ca n A 8927 lines compared with the observed spectrum of 
the Kurucz et al. (1984) solar flux atlas (bold dots). The fitting 
parameters are presented in Table@] It should be stressed that 
different Ca abundances are required to achieve the best fits 
of the observed profile in the NLTE (Sh — 0) and LTE cases. 
Everywhere V m ; c = 0.9kms _1 
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Fig. 6. Solar NLTE (Sn = 0) Ca abundances from the Cai 
(filled circles) and Ca n (open diamonds) lines plotted as func- 
tion of the line equivalent width. Top, middle, and bottom pan- 
els show results based on the measured fy, OP, and NIST data, 
respectively. In each panel, the mean abundance derived from 
the Ca i lines is shown by the dotted line and the mean values 
and their standard deviations are quoted for Ca i and Ca n 



A 8927 of Ca n. It is clearly seen that the observed line core of 
Ca i A 6439 is asymmetric. This is, most probably, due to atmo- 
spheric inhomogeneity and, therefore, cannot be reproduced in 
the framework of a ID analysis. 

Using the obtained values log g/s2 and fy from different 
sets, we compute Ca abundances from individual lines. Fig.[6] 
illustrates results for NLTE (S H =0) calculations. It can be 
seen that the laboratory oscillator strengths provide the high- 
est accuracy in the absolute solar abundance derived from the 
lines of Ca i. The NLTE and LTE averages from 23 lines are 
logs C ai(LAB, S n = 0) = 6.37 + 0.06, logs C ai(LAB, S H = 
0.1) = 6.36 + 0.06, loge C ai(LAB, S H = 1) = 6.35 + 0.06, and 



logecai(LAB, LTE) = 6.34 + 0.06. Throughout this paper, stan- 
dard deviations are quoted. NLTE effects for Ca i in the Sun are 
very small, and further we present the results only for Sn = 0. 
Ca abundances based on OP fy for 20 Ca i lines reveal a large 
spread of data with the mean values log ecai(OP, Sn — 0) = 6.49 
± 0.23 and loge C ai(OP, LTE) = 6.46 + 0.21. Excluding the one 
line, A 5867, with the largest contribution to standard deviation, 
we obtain loge Ca i(OP, S H = 0) = 6.45 ± 0.14 and loge Ca i(OP, 
LTE) = 6.43 + 0.14. From 20 Cai lines with the NIST data 
available, the averages equal log ec a i(NIST, Sn = 0) = 6.29 + 
0.13 and log e Ca i(NIST, LTE) = 6.26 ± 0.13. The results based 
on the VALD data are close to that for laboratory measure- 
ments because, for most lines investigated, VALD's oscillator 
strengths have been taken from Smith & Raggett ( 1981 1. Tests 
show that the NLTE abundance from Ca i lines is not sensitive 
to a variation of electronic collisional rates in SE calculations. 
Using the formula of van Regemorter ( 1962 1 and IPM (Seaton 
1 1962 a I data leads to the mean Ca abundances, consistent within 
0.01 dex. In each case, either NLTE or LTE, for any set of fy, 
no significant correlation of individual Ca abundance is found 
with the line strength. 

Only predicted oscillator strengths are available for the Ca n 
lines which are used to determine the absolute solar abundance 
from Can lines. OP data provide high accuracy for the desired 
value provided that an SE approach is applied. NLTE correc- 
tions remove a trend of the Ca n abundance with line strength 
obtained under the LTE assumption. Averages from eight lines 
are loge Ca ii(OP, S H = 0) = 6.38 + 0.07, log£ C aii(OP, S H = 
0.1) = 6.40 ± 0.06, and loge C aii(OP, S H = 1) = 6.43 + 0.08. It 
should be noted that a variation in log C& values within 0.6 dex 
for the lines of multiplets 5p — 5d and Ad - 4/ sensitive to 
van der Waals broadening changes the mean Can abundance 
by only 0.01 dex. The NIST database contains only three Can 
lines of our interest and their fy are very close to the corre- 
sponding values from OP calculations. For the VALD data, the 
NLTE and LTE mean Ca n abundances equal log ecaii( VALD, 
S H = 0) = 6.65 + 0.32 and log e C aii( VALD, LTE) = 6.74 + 0.29. 
The large standard deviation is mainly caused by A 6456 and 
A5339. In contrast to Cai, NLTE Can abundance is sensitive 
to a variation of electronic collisional rates in SE calculations. 
The value determined using the formula of van Regemorter 
(1962} and assuming Sn — is larger by 0.05 dex compared 
to the corresponding value obtained for the standard recipe of 
electronic collisional rates. 

Taking into account the highest accuracy of the mean val- 
ues, we prefer to use as final values the Ca n abundance based 
on OP oscillator strengths and Ca i abundance based on the lab- 
oratory fy. They agree within < 0.04 dex provided that Sn 
< 0.1 and the standard recipe of electronic collisional rates 
is applied. In any other case, the difference equals 0.08 dex 
(Sh - 1 or van Regemorter's electronic collisional rates) to 
0.12 dex (LTE). Thus, using the theoretical MAFAGS model 
atmosphere, we find solar Ca abundance to lie between 6.36 
and 6.40. 

We obtain here also the fitting parameters of solar Cai 
A 4226 and Ca n A 8498 lines necessary for further analysis of 
metal-poor stars. For both lines, oscillator strengths are taken 
from laboratory measurements (Smith & Gallagher 1966] for 
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Table 4. Atomic data for the Ca i and Ca n lines: column 2 contains the multiplet numbers accordingly to Moore (1972); columns 
4-7 give log gf from various sources including laboratory measurements (column LAB); log C(, of individual lines are from the 
sources cited in the bottom part of this table except for the unmarked values which are based on the Smith (1981) data. Solar 
l°gg/ e Ca values were determined from NLTE analysis of the Ca line profiles in the Kitt Peak Solar Atlas (Kurucz et al. 1984) 
neglecting hydrogenic collisions (Su = 0). Columns 10 - 12 present the abundance corrections A x = log ex - loges H =o where 
X takes the meaning 0.1, 1, and LTE for Sh = 0.1,Sh = 1 an d the LTE assumption, respectively. The last column contains solar 
equivalent widths 
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A 4226; Theodosiou [T9891 for ,18498), and the van der Waals 
damping constants are based on A&O'M's, data (Table^J. For 
A 8498, isotope structure is taken into account with atomic data 
from Tabled A good fit of the observed solar flux profile of 
/14226 in the wings is achieved for log eca = 6.29 when hydro- 



genic collisions are ignored and for logeca = 6.21 in the LTE 
case. NLTE effects are negligible for the A 8498 line wings. The 
best fit is achieved at a Ca abundance of log eca = 6.28. 

Recent determinations of solar photospheric Ca abundance 
based on a 3D LTE analysis (Asplund et al. |2005 ) give "ex- 
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cellent agreement between the two ionization stages" and the 
average of the two: logeca = 6.31+0.04. The average ob- 
tained from our ID NLTE analysis is log eca = 6.38 ± 0.06. We 
showed above that the absolute Ca abundance depends on the 
adopted values of oscillator strengths and van der Waals damp- 
ing constants. The cited paper of Asplund et al. does not give 
enough information for a discussion of possible sources of the 
found discrepancy. 

5. Stellar sample, observations and stellar 
parameters 

Our sample consists of eight stars. HD 61421 (Procyon) is se- 
lected as a fundamental star with nearly the same list of de- 
tected Ca lines as in the Sun. It is especially important for 
checking the Can high excitation lines. The other objects are 
metal-poor stars with metallicity ranging between [Fe/H] = 
-1.35 and [Fe/H] = -2.43. They give an opportunity to study 
a formation of some of the strongest Ca lines, Ca i A 4226 and 
Ca ii A 8498 which presumably become of purely photospheric 
origin in this metallicity range. 

Four metal-poor stars, HD 29907, HD 59392, HD 140283 
and BD -4°3208 were observed using the Ultraviolet and 
Visual Echelle Spectrograph UVES at the 8m ESO VLT UT2 
telescope on Cerro Paranal. At least two exposures were ob- 
tained for each star. Spectral resolving power is about 60 000. 
The data cover an approximate spectral range of 4000 - 7000 A. 
The signal-to-noise ratio is 200 or higher over the whole spec- 
tral range. No near-IR spectra for HD 29907 and HD 59392 are 
available to us. 

Observational data for all other stars (and for Ca n 8498 in 
BD -4°3208) are taken from Korn et al. 120031 with overall 
similar data-quality specifications as for the UVES data. For 
these data, the wavelength coverage is 4200 -9000 A. High- 
quality observations for Can 8498 in HD 140283 were taken 



from the ESO UVESPOP survey (Bagnulo et al. 2005 ). 

We use stellar parameters determined in our earlier stud- 



ies (Korn et al. f2003i Mashonkina et al. 2003). In short, ef- 
fective temperatures were determined from Balmer line pro- 
file fitting, Ha and H/3, with a statistical error estimated at the 
level of 70 K, surface gravities from the Hipparcos parallaxes 
with masses determined from the tracks of VandenBerg et al. 
(20001. The errors obtained from adding the squared errors of 
parallax and mass are quoted in Table[5] The iron abundance 
and microturbulence velocity were obtained requiring the de- 
rived [Fe/H]n 2 abundances not to depend on line strength. For 
metal-poor stars, a-enhanced models are calculated with abun- 
dances of a-elements O, Mg, Si and Ca scaled by the stellar 
Mg/Fe ratio. The latter is determined in this study from NLTE 
analysis of four Mg i lines, A 457 1 , A 4703, A 5528, and A 57 1 1 . 
We use the model atom and atomic data for Mg i described by 
Gehren et al. (2004). Stellar parameters are given in Table|5] 

The microturbulence values are somewhat lower than val- 
ues previously published by us as we now utilize the broad- 



2 The subscript in [X/Y]n indicates that the abundance for element 
X is determined from the singly ionized species. Likewise, a subscript 
I refers to neutral lines used in the abundance determination. 
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Fig. 7. Trends of abundance with line strength determined from 
the Ca i lines in Procyon for our best NLTE model (top) and un- 
der the assumption of LTE (bottom). Note the steep trend of the 
LTE [Ca/H] values with line strength above Wa - 100 mA. The 
mean values given are computed from weak lines only (W^ < 
100 mA). 



ening parameters of Barklem & Aspelund-Johansson (2005) 
which play a role even in the metal-poor stars as the analysis 
is differential to the Sun. In all metal-poor stars, the discrimi- 
nation of the microturbulence values relies heavily on the three 
strong Fen lines of multiplet 42 (AA 4923, 5018 and 5169 A). 
All three lines seem to require substantially higher microtur- 
bulence values in the Sun. Enforcing a fixed solar microturbu- 
lence of y m ; c = 0.9 km s -1 thus means that they cannot be well 
fitted. We estimate that, depending on the choices made in the 
fitting of solar lines and the line selection for the metal-poor 
stars, microturbulence values can easily vary by 0.2 km s -1 . In 
particular, lower values seem possible (e.g. disregarding Fe n 
5169 A). This is important to bear in mind for the discussion of 
the strong lines in Sections[6]and0. 

Our results are based on line profile analysis. In order to 
compare with observations, computed synthetic profiles are 
convolved with a profile that combines instrumental broaden- 
ing with a Gaussian profile and broadening by macroturbulence 
with a radial-tangential profile. Only slow rotators are included 
in our sample, Vsini < 1.5kms _1 except for Procyon with 
V sin/ = 2.6kms~ 1 (Fuhrmann 1998) and therefore the rota- 
tional velocity and macroturbulence value cannot be separated 
at the spectral resolving power of our spectra. We thus treat 
their overal effect as radial-tangential macroturbulence. Only 
for Procyon, rotational broadening and broadening by macro- 
turbulence are treated separately. The macroturbulence value 
was determined for each star in our previous studies from the 
analysis of an extended list of lines of Fei/n, Mgi, etc. Here, 
Vmac was allowed to vary by +0.3 km s (l<x). 

6. The Ca i versus Ca n in the selected stars 

In this section, we test NLTE formation of the Can lines in 
two steps. In the first one, the Can high excitation lines are 
examined in the solar metallicity star Procyon. These lines are 
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Table 5. Stellar parameters and obtained NLTE elemental abundance ratios of the selected stars. V m i C is given in km s _I . 
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more sensitive to details of NLTE calculations than A 8498 (see 
Table[3}. We then derive element abundances from the Ca i lines 
and from Ca n A 8498 for five metal-poor stars and inspect the 
difference Aloge(Cai - Can). In the temperature regime we 
are concerned with and at metallicities [Fe/H] between -2 and 
-3, the equivalent width ratio W(Ca i)/W(Ca n 8498) is weakly 
sensitive to a variation of T e g and log g, independent of what 
Ca i line is taken. For example, log W(Ca i 6439) / W(Ca n 
8498) = -0.88, -0.81, -0.77, and -0.78 for the models with 
Teff I logg / [Fe/H] = 5000K / 3.0 / -2, 5000K / 4.0 / -2, 
6000 K / 3.0 / -2, and 6000 K / 4.0 / -2, respectively. This 
can be understood because the van der Waals broadened line 
wings give a significant contribution to W(Ca n 8498) even at 
[Ca/H] = -2.6 (for the models with [Fe/H] = -3). When logg 
increases, strengthening of the A 8498 line wings is compen- 
sated by weakening NLTE effects. In such conditions, a dif- 
ference between Ca abundances determined from two ioniza- 
tion stages, if present, will point to shortcomings of the NLTE 
treatment of Ca lines rather than to the uncertainty of stellar 
parameters. 

Analysis of stellar spectra is made line-by-line differen- 
tially with respect to the Sun. At the LTE assumption, the cores 
of many lines cannot be fitted. In such cases, Ca abundance is 
derived from the line wing fitting. 

Procyon 

Irrespective of the assumed efficiency of hydrogen collisions, 
average NLTE effects are very small among weak Cai lines 
(W A < 100mA) and vary between -0.02 dex (S H = 1) and 
-0.03 dex {Su — 0). However, the SE approach is able to re- 
move a steep trend with line strength among strong Ca i lines 
seen in LTE (see Fig.0. For individual strong lines (like Cai 
6439 A) NLTE corrections exceed -0.3 dex. Among weak lines 
(Wa < 100 mA), the line-to-line scatter is reduced by all three 
NLTE models and attains its minimal value (<x = 0.027 dex) 
at S h = 0.1. Surprisingly, the mean abundance derived from 
Cai lines is clearly subsolar, [Ca/Hh = -0.12 + 0.04 (LTE) to 
-0.15 + 0.04 (S H = 0). 

Abundances for Ca n are derived from six high-excitation 
lines. The mean NLTE value is [Ca/H]„(NLTE) = 0.07 + 0.05 
independent of the assumed strength of hydrogenic colli- 
sions. LTE abundances show a mean value of [Ca/H]n(LTE) = 
0.14 + 0.05. An application of the formula of van Regemorter 
( 1962 1 for calculation of electronic collisional rates results in 
an even larger difference A log e(Ca i - Ca n). 




8496 8497 8498 8499 8500 



Fig. 8. The best NLTE (S H = 0.1) fits (continuous line) of the 
observed profiles of Can A 8498 in spectra of HD 19445 and 
HD 103095 (bold dots). The obtained NLTE Ca abundances are 
presented in Table [5] For HD 19445, line profiles are shifted 
along Y axis 

Can the stellar parameters be the cause of this discrepancy? 
Allende Prieto et al. (2002) estimate an error of the effec- 
tive temperature of 49 K. Ramirez & Melendez ( 2005 1 revise 
the fundamental effective temperature of Procyon to rf™ d = 
6532 + 39 K and derive 7^™ = 6591 ± 73 K using their in- 
frared flux method. Assuming T e ff = 6590 K (a value favoured 
by Korn et al. 2003 for the NLTE ionization equilibrium of 
Fe i/ii) and keeping the values of log g, [Fe/H] and V m j c un- 
changed we obtain mean values [Ca/H]i(NLTE) = -0.10 ± 0.04 
and [Ca/H] n (NLTE) = 0.05 + 0.05. As the discrepancy lies with 
weak lines, modifications to the microturbulence value cannot 
remove this discrepancy, either. Aufdenberg et al. ( 2005 ) con- 
clude that the interferometric data provide evidence for con- 
vective overshooting in Procyon's atmosphere. We checked a 
model which includes an "approximate overshooting" prescrip- 
tion for convective flux transport in a mixing-length formal- 
ism according to Castelli et al. ( 1997). The investigated lines 
of both Ca i and Ca n become weaker in this model compared 
to our standard MAFAGS model, and the derived abundances 
of Cai and Can increase by 0.08 - 0.10 dex for different lines. 
However, the difference Aloge(Cai - Can) remains at the 
same level. 

The effect of atmospheric temperature inhomogeneities on 
Ca i/ii is expected to be of the same order of magnitude as that 
for neutral and singly ionized iron lines calculated by Allende 
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Prieto et al. (2002i. They show that weak lines (W A < 50 mA) 
of both ionization stages, Fe i and Fe n, are weakened com- 
pared to a classical ID analysis, such that the derived Fe abun- 
dance increases by 0.05 dex and 0.04 dex, respectively. Three- 
dimensional simulations might lead to the higher Ca abundance 
for Procyon compared to our results, however, the discrepancy 
between Ca i and Ca n is unlikely to be removed. Thus, while 
our NLTE analysis leads to a better agreement of the Ca abun- 
dance from the two ionization stages, the discrepancy exceeds 
3cr under standard assumptions about Procyon's stellar param- 
eters. 
HP 19445 

The abundance is determined from 15 Cai lines with equiva- 
lent widths between 4 mA and 63 mA. The accuracy of the ob- 
tained differential abundances is at the level of cr = 0.04 dex for 
Sh = 0, 0.1 and the LTE case and slightly worse, <x = 0.06 dex 
for S h = 1 ■ The difference between NLTE and LTE Ca i abun- 
dances depends strongly on the assumed value for S h- It equals 
+0. 13 dex for Sh — and reduces down to +0.03 dex for Sh = 
1 . In contrast, NLTE abundances from Ca n A 8498 are smaller 
than in LTE. The best fit of A 8498 achieved at S H = 0.1 and 
[Ca/Fe] = 0.40 is shown in Fig. [3] The difference between 
Ca i and Ca n abundances is found to be 0.04 dex, -0.04 dex, 
-0.11 dex and -0.19 dex for S H = 0, S B = 0.1, S H = 1 andLTE, 
respectively. An uncertainty in microturbulence of 0.2kms _1 
hardly affects the Ca i abundance and leads to a change in Ca n 
of 0.03 dex. This does not destroy the agreement of the NLTE 
Ca i and Ca n abundances if S h = or 0. 1 is assumed. 
HP 84937 

The abundance obtained from Ca i lines vary between [Ca/Fe]i 
= 0.47 + 0.04 and [Ca/Fe]i = 0.39 + 0.04, when changing 5 H 
between and oo (LTE). Similarly to HP 19445, the abun- 
dances from the two ionization stages agree much better in the 
NLTE case than in LTE one. Alogs(Cai - Can) = +0.11 dex 
and -0.02 dex for NLTE abundances at Sh =0 and 1, re- 
spectively, while the difference of LTE abundances equals - 
0.23 dex. The best agreement is found for an Sh value between 
0.1 and 1. 
HP 103095 

The determined NLTE abundances range between [Ca/Fe]i = 
0.37 + 0.04 and [Ca/Fe]: = 0.29 + 0.04, depending on the as- 
sumed Sh value. NLTE effects for Can /18498 occur only in 
the very core, and the same Ca abundance is obtained with 
[Ca/Fe]n = 0.29, independent of the used theory of line forma- 
tion. The best fit of this line is achieved for Sh =0.1 (Fig.|8j. 
The difference between Ca i and Ca n is within the mutual error 
bars in all cases: Aloge(Cai - Can) = +0.08 dex, +0.01 dex, 
0.00 dex, and +0.03 dex for S H = 0, 0.1, 1, and LTE, respec- 
tively. It should be noted that the discrepancy is larger when 
hydrogenic collisions are neglected (cf. HP 84937 above). 
HP 140283 

The [Ca/Fe]: ratios vary between +0.16 + 0.03 (LTE) and 
+0.33 + 0.03 (Sh — 0) and are determined from six lines with 
equivalent widths between 6 and 40 mA. At a microturbulence 
value of V m ic = 1.65 km s~' the Cai resonance line (W,\ = 146 
mA) is in good agreement with the weak lines under the as- 
sumption of LTE, but yields too low abundances for all NLTE 
models. We comment on this behavior in the next section. Like 



in the case of HP 84937, A log e(Ca i - Ca n) vanishes between 
Sh = 0.1 and 1. In LTE, however, its value is -0.3 dex. 
BP -4°3208 

This star is found to be very similar to HP 84937, both with 
respect to its stellar parameters and its behavior in terms 
of Aloge(Cai - Can). The abundances obtained from 14 
Cai lines with equivalent widths between 4mA and 40mA 
give mean values of [Ca/Fe]i(5 H = 0.1) = 0.43 + 0.04 and 
[Ca/Fe]:(LTE) = 0.34 + 0.04. Good agreement between Cai 
and Can is again found for an Sh range of to 0.1, while the 
difference increases above Sh = 0.1 and reaches -0.35 dex in 
the case of LTE. 

Summarizing our results for the five metal-poor stars, we 
conclude that, within the modelling uncertainties, NLTE leads 
to consistent Ca abundances derived from the two ionization 
stages while LTE fails to give consistent results. We find that 
ignoring hydrogenic collisions results in too strong NLTE ef- 
fects for Cai/n. Based on the results obtained for the Sun and 
the stars presented in this section, our best choice for a scaling 
factor applied to Steenbock & Holweger's ( 1984 > formula is 
Sh =0.1. Final Cai and Can abundances corresponding to Sh 
= 0.1 are presented in Table[5] 

7. The Ca i resonance line versus subordinate 
lines in metal-poor stars 

Cai A 4226 is the only neutral calcium line which can be de- 
tected in extremely metal-poor stars. In this section, we test 
NLTE formation of this line using the metal-poor stars where, 
on one hand, the resonance line becomes, probably, of purely 
photospheric origin and, on other hand, Ca i subordinate lines 
are still measurable and provide a reliable value of the Ca abun- 
dance. In addition to the five metal-poor stars discussed above, 
another two stars are studied here. We start by determining the 
[Ca/Fe]: ratio for them. 
HP 29907 

The Cai abundance is determined from 16 lines with equiv- 
alent widths between 18 mA and 174 mA. NLTE removes a 
trend with line strength displayed by LTE abundances and re- 
sults in the smaller standard deviation, cr = 0.05 dex, compared 
to cr = 0.08 dex in the LTE case. Mean NLTE (S H =0.1) and 
LTE abundances agree within 0.01 dex. 
HP 59392 

The NLTE and LTE mean abundances derived from 17 
Cai lines turn out quite similar with [Ca/Fe] (S h = 0.1) = 
0.26 + 0.04 and [Ca/Fe] (LTE) = 0.24 + 0.05. 

Our calculations show that, for every star, overionization 
in the atmospheric layers below logrsooo = -1 (Fig.|2j leads 
to weakening the line wings of A 4226 compared to LTE 
case, while a steep decrease of the departure coefficient ratio 
b{4p)/b(4s) above that depth point results in opposite effect 
for the line core. In the two coolest stars of our small sam- 
ple, HP 29907 and HP 103095, the Ca i resonance line is very 
strong with the core formed in the uppermost atmospheric lay- 
ers near logrsooo = -4.4 and -4.9, respectively (we indicate 
here a location of line center optical depth unity). However, 
total line absorption is dominated by the van der Waals broad- 
ened wings and not by the core. For both stars, we find that the 
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Fig. 9. Theoretical NLTE (Sh — 0.1, continuous line) profiles 
of the Ca i A 4226 line compared with the observed spectra of 
HD 29907 and HD 103095 (bold dots). For each star, the Ca 
abundance adopted in the NLTE calculations is taken from the 
analysis of the Cai subordinate lines (see Table [5}, and this 
abundances allows to fit the resonance line as well. Offsets in y 
are applied for better illustration 

NLTE theoretical profile computed with Ca abundance deter- 
mined from the Cai subordinate lines describes the observed 
profile well except for the very core which can be influenced 
by the star's chromosphere. Results are illustrated in Fig. [9] 

For HD 84937 and HD 19445, the agreement between sub- 
ordinate lines of Cai, Can 8498 and Cai 4226 is excellent 
when Sh- 0.1 is assuming. Cai 4226 requires a [Ca/Fe] ratio 
of 0.41 for HD 84937 and of 0.36 for HD 19445. 

A less consistent picture emerges for the remaining three 
stars. When Ca abundance is fixed at the value derived from 
the analysis of the Ca i subordinate lines, the half- width of the 
theoretical NLTE profile of A 4226 is obtained to be larger than 
the observed one. To fit the observed line width, a smaller Ca 
abundance is required: by 0.08 dex for HD 59292, 0.09 dex for 
HD 140283, and 0.22 dex for BD-4°3208. Discrepancies are 
smaller in LTE for these stars, but larger for e.g. HD 84937. 

In each star investigated, Cai 4226 lies on the saturated 
part of the curve of growth, and such discrepancy could 
well be related to uncertainties in the microturbulence value: 
AVmic = -O.lkms -1 translates to A log e(A 4226) = +0.03 dex 
to +0.04 dex for different stars. However, for a star like 
BD-4°3208, a substantial reduction of the microturbulence 
would be required, which in turn would affect the good agree- 
ment between Ca i and Ca n reported above. We note in passing 
that a similar discrepancy is found for subordinate lines of Mg i 
and the Mgib triplet lines in these stars. 

Are these then shortcomings of our SE calculations? 
Collisions are inefficient in statistical equilibrium of atoms in 
the layers, where the Doppler core of A 4226 is formed, be- 
tween logT5ooo = -1 and -3, for the most metal-poor stars of 
our sample, BD -4°3208 and HD 140283. The related non-LTE 
effect is therefore entirely due to photoionization which should 
be modelled accurately (given the correctness of the OP pho- 
toionization cross-sections.) The explanation can lie with the 
adopted one-dimensional atmospheric models. Based on the re- 
cent results of Shchukina et al. ( 2005 ) for weak and moderately 



strong (W A < 80 mA) lines of Fei in HD 140283, we expect 
an overall small effect of atmospheric-temperature and veloc- 
ity inhomogeneities on the Ca i subordinate lines in our sample 
of metal-poor stars. The Cai resonance line may be a differ- 
ent case because it is formed over the more extended range of 
atmospheric depths. For a fully quantitative understanding of 
its formation in cool stars, further studies are required on the 
basis of advanced model atmospheres. One clear advantage of 
such studies will be the removal of adjustable parameters like 
microturbulence from which the current modelling potentially 
suffers. 

8. The Ca i/Ca n ionization equilibrium as indicator 
of surface gravity for extremely metal-poor stars 

In this section, we consider the possibility of using the W(Ca i 
4226) / W(Can 3933) and W(Cai 4226) / W(Can 8498) equiv- 
alent width ratios as indicators of surface gravity for extremely 
metal-poor stars. At extremely low Ca abundance, the line 
wings do not contribute anymore to W(Can 8498). The line 
is expected to be strengthened with decreasing logg due to 
decreasing H continuous absorption. NLTE effects for Can 
A 8498 are amplified in the same direction. Thus, the W(Ca i 
4226) / W(Can 8498) ratio should increase when logg goes 
up. The Ca n resonance lines remain strong even at [Ca/H] = 
-5, and their van der Waals broadened wings are weakened 
with decreasing log g. The Ca i resonance line is weakened in 
the same direction due to amplified overionization. Thus, the 
W(Cai 4226) / W(Can 3933) ratio is rather insensitive to vari- 
ation of log g. Both ratios were calculated for a small grid of 
models with [Fe/H] = -4.34; the Ca abundance was adopted to 
be [Ca/H] = -4.9. The results corresponding to T e g = 5500/f 
are plotted in Fig. EH as a function of logg. It is clear that the 
W(Ca i 4226) / W(Ca n 8498) ratio can be used to determine 
the surface gravity for extremely metal-poor stars, contrary to 
the ratio involving the Ca n resonance line(s) which is nearly 
constant over the logg range between 2.5 and 4.5. An applica- 
tion of this technique to the two known ultra-metal-poor stars 
(HE 0107- 5240, Christlieb et al.|2002 HE 1327-2326, Frebel 
et al. 2005 1 will be presented in a forthcoming paper. 

9. Concluding remarks 

In this study, NLTE line formation of an extended list of Ca i 
and Can lines was considered for the temperatures ranging 
between T e ff = 5000/T and T e ff = 6000^, surface gravities 
logg = 3.0 and 4.0, and metallicities from [Fe/H] = down to 
[Fe/H] = -4.34. For every Ca i line, departures from LTE affect 
significantly its profile over the whole range of stellar param- 
eters. If a line is strong (multiplets 2, 3, 4, 18, 21, and 47 at 
[Fe/H] > -2 and the remaining multiplets at [Fe/H] > — 1), its 
wings are weakened but the core is strengthened compared to 
the LTE case. The value and sign of NLTE abundance correc- 
tion are defined by a relative contribution of the core and the 
wings to the overall line strength. When the line becomes weak 
due to decreasing Ca abundance, NLTE leads to depleted total 
absorption in the line and positive abundance correction. NLTE 
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Fig. 10. The Cai 4226 / Can 3933 and Cai 4226 / Can 8498 
equivalent width ratios as a function of the surface gravity for 
the models with r eff = 5500 K and [Fe/H] = -4.34. A Ca abun- 
dance of [Ca/H] = -4.9 is adopted everywhere 

effects are very large at extremely low Ca abundance. For ex- 
ample, at [Ca/H] = -4.9 NLTE abundance correction for Cai 
A 4226 can exceed 0.5 dex. Thus, for a given line, Anlte de- 
pends on r e ff, logg, [Ca/H], and microturbulence value, and, 
for a given model, Anlte is different in value and sign for a va- 
riety of Cai lines. Any interpolation of NLTE results for Cai 
has to be performed with caution. A different situation is found 
for Ca ii. For every Ca n line, NLTE leads to enhanced absorp- 
tion in the line core and negative abundance correction over 
the whole range of stellar parameters. The absolute value of 
Anlte is defined by the relative contribution of the core and the 
wings to equivalent width. For example, NLTE corrections re- 
main very small, < 0.02 dex, for the resonance lines and grow 
in absolute value with decreasing Ca abundance for the IR lines 
of multiplet 3d - 4p exceeding 0.4 dex in a metal-poor model 
with r eff = 6000K, logg = 3.0 and [Fe/H] = -3. Such correc- 
tions are important to be considered for the correct interpreta- 
tion of the near-IR spectroscopic data collected in the RAVE 
survey and the future ESA Gaia satellite mission. 

Empirical evidence is found from the analysis of stellar 
spectra that inelastic collisions with hydrogen atoms serve as an 
additional source of thermalization in cool (metal-poor) stars. 
Our best choice for a scaling factor to the formula of Steenbock 
& Holweger ( 1984) for calculation of hydrogenic collisions is 
Sh = 0.1. Disregarding hydrogenic collisions completely will 
lead to overestimated NLTE effect at low metallicity. 

Taking advantage of our SE approach and accurate atomic 
data for the investigated lines, we obtain good agreement, 
within 0.04 dex, between absolute solar Ca i (from 23 lines) and 
Ca ii (from 8 lines) abundances, and the average of the two is 
loge° a = 6.38 + 0.06. 

Calcium abundances from two ionization stages are exam- 
ined for the first time for five metal-poor stars. We show that 
NLTE largely removes obvious discrepancies between Ca i and 
Ca ii obtained under the LTE assumption. 

Solving the restricted NLTE problem for calcium in "clas- 
sical" one-dimensional LTE model atmospheres, we qualita- 
tively understand the formation of the Cai resonance line in 



metal-poor stars, where it is of purely photospheric origin. In 
agreement with observations, NLTE predicts weakening of the 
line wings and strengthening of the line core compared to the 
LTE case. Residual discrepancies may be related to the use of 
classical model atmospheres with adjustable parameters like 
microturbulence . 
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